Extended Bose-Hubbard model with incompressible states at fractional numbers 
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The Bose-Hubbard model is extended to include nearest and far neighbor interactions and is 
related to the fractional quantum Hall effect (FQHE). Both models may be studied in optical 
lattices with quantum gases. The ground state is calculated for the extended Bose-Hubbard model 
with strong repulsive interactions (weak hopping). Incompressible Mott insulator states are found 
at rational filling fractions compatible with the principal and secondary FQHE filling fractions of 
the lowest Landau levels observed experimentally. It is discussed to which extent these states at 
fractional filling survive or undergoes a Mott insulator transition to a superfluid as hopping terms 
are included. 

PACS numbers: 03.75.Hh73.43.-f,71.10.-w,03.75.Lm 
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I. INTRODUCTION 

Quantum gases provide new insight in strongly corre- 
lated quantum systems and phase transitions in systems 
of interacting bosons and fermions. A variety of interest- 
ing phases exist in various dimensions such as superflu- 
ids, Mott insulators, charge density waves, incompress- 
ible states in the integer and fractional quantum Hall ef- 
fect (FQHE) U0HS@i,etc. Lattices appear naturally 
in solids and spin models laSSQi an d recently also op- 
tical lattices with Bose atoms have been studied 0, ^3 ■ 
The su perf luid to Mott insulator transition has been ob- 
served |Ljj as predicted from the Bose-Hubbard model. 
It has been suggested that by applying an additional 
magnetic field the optical lattices provide a clean system 
without impurities, where the FQHE can be studied in 
detail. Alternatively, rapidly rotated harmonic oscillator 
traps with cold Fermi or Bose atoms may also form a 2D 
Hall fluid [H E E3 • 

These possibilities stimulate renewed interest in the 
long standing problem of the FQHE and its numerous in- 
compressible quantum liquid states at mysterious filling 
fractions of the lowest Landau levels 0,0- ln a recent ex- 
periment Q new FQHE fractions are found at 3/8, 4/11, 
5/13, 6/17, 7/11, etc., which fall outside the Laughlin 
fractions l/(2m + 1) and Jain fractions m/(2m ± 1) for 
integer m. 

The Hubbard model is widely used to describe a num- 
ber of systems such as strongly correlated electrons. 
Quantum gases of Bose atoms in optical lattices have re- 
cently been described by the Bose-Hubbard model, which 
in one dimension reads 
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BH 
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Here, t\ is the hopping parameter to nearest site and 
Uq is the on-site interaction energy. The system is a su- 
perfluid except in the cases where the filling fraction is 



precisely an integer and the hopping term is sufficiently 
small, where a transition occurs from a superfluid to a 
Mott insulator as observed in |l3j|. Including nearest 
neighbor interaction leads to charge density waves and 
crystalline structure at half filling and a related Mott 
transition [l7l ITS Il9j . 

We shall be interested in extended Bose-Hubbard mod- 
els (EBHM) including not only on-site interactions Uq 
but also interactions U n between particles n — 1, 2, 3, 4, ... 
sites away. In the strongly interacting (atomic) limit or 
equivalently weak hopping, t\ <C U„ , the extended Bose 
and Fermi Hubbard models reduce to well known spin 
models. The Ising model [6| has only nearest neighbor 
U\ interactions but next-nearest neighbor interactions U<i 
have also been included (see, e.g., H,@)- 

The purpose of this work is to point to a connection 
between the Hamiltonian for the FQHE and the EBHM 
as outlined in section II. Subsequently in section III to 
find the ground state of the EBHM in the strongly in- 
teracting limit, calculate its energy, chemical potential, 
pressure, etc. Most importantly, a spectrum of incom- 
pressible Mott insulator states are found whenever the 
filling fraction is a rational number q/n. In section IV 
the Bose filling fractions are related to the fermion FQHE 
case, and it shown that they contain the Laughlin and 
Jain sequences 0, |2(j as well as all the secondary FQHE 
fractions found experimentally |3j. Finally in section V 
the effect of hopping terms and the survival of Mott in- 
sulator states at fractional filling is discussed. 



II. APPROXIMATING THE FQHE AS AN 
EXTENDED BOSE-HUBBARD MODEL 

The fractional quantum Hall systems [j| can be es- 
sentially reduced to the problem of N spinless elec- 
trons with repulsive interactions in M available states 
[22|. In the disk geometry the states are ipk(z) — 
z k exp(-|z| 2 /A)/V2TT2 k kl, k = 0, 1, M - 1; 
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iy = re 1 ^ in units of the magnetic length I |2]| . 

When the filling fraction is exactly v — N/M = 1, 
1/3, 1/5, 1/7, etc., the ground state is the antisymmetric 



2 



Laughlin wave function 



N 



N 



(2) 



with m = 1, 3, 5, 7, etc., and normalization factor A 2 . 
Analogously, a symmetric Laughlin wave function with 
m = 2,4,6,... can be separated for bosons instead of 
fcrmions. 

In the lowest Landau level (N < M) the total wave 
function also contains a symmetric part ^ s 



* = * S (zi,...,Zjv)*r, 



(3) 



The remaining FQHE problem is to find the incompress- 
ible states of VPs besides the fractions v = 1, 1/3, 1/5, ... 
implicit in the Laughlin wave function. At these fractions 

*5 - 1. 

In the following we shall concentrate on the lowest Lan- 
dau level, N < M , and study the symmetric wave func- 
tion ^5 for N particles in Mb — M — mN available 
states. Here the number of states, that the Laughlin 
wave function use up, have been subtracted. By this 
"bosonization" the N particles has thus become N spin- 
less bosons in Mb states with a symmetric many-boson 
wave function ^ $ (see, e.g. |23| for details on bosoniza- 
tion). The Bose filling fraction vb = N/Mb is related to 
the fermion one by 



N 
M 



vb 



1 + mvB 



(4) 



with m = 1, 3, 5,.... 

We can choose any other orthogonal basis from tpk(z). 
This may be convenient if a many-body wave function 
can be constructed from this new basis such that the 
Hamiltonian approximately diagonalizes. This would 
greatly simplify the search for the ground state of the 
FQHE. A candidate for such a convenient basis are the 
localized Wannier-type wave functions 



1 



y/2¥Ml 



Mb-1 

£ Mz)e~ i2 * kj/MB ■ (5) 



fe=0 



where j — 0,1,..., Mb — 1. As described in Appendix 
A these wave functions are localized in angle <p around 
2irj/M B . 

Dyakonov argues in Ref . |22j the similar localized wave 
functions are useful for constructing a crystal-like many- 
body state that contains at least some of the right fea- 
tures of the FQHE ground state. The localized wave 
functions of Ref. (2^ are an approximate model for the 
FQHE on a circle. It is simple and illustrative, and is 
therefore described in Appendix A. 

It should be emphasized that the localized states are 
simply another basis. They are not generated by an 
external periodic potential as, e.g., an atomic lattice. 
Whether the basis is "better" will depend on whether 



the repulsion between particles dominates over hopping 
leading to localized vs. conduction states respectively. It 
will be argued below that localized states of Eq. JSJ may 
be the preferred ones. 

The Hamiltonian in the lowest Landau level is now a 
one dimensional lattice model containing only repulsive 
two-body interaction V between the particles 



H 



E 

l<i<j<N 



V(zi - z s ) . 



(6) 



The interactions can in the following be almost any re- 
pulsive potential including the standard Coulomb: V — 
e 2 l\z\ — z%\. The delta function potential should be 
treated with care as discussed in Appendix A. 

The kinetic energies and other degrees of freedom are 
absorbed in the cyclotron energies, huj c = eB/mc, of the 
lowest Landau level. The cyclotron energies have been 
subtracted so that only interaction energies remain in 
the Hamiltonian. In this sense the system is strongly 
interacting. All energies of this Hamiltonian will scale 
with the overall strength of the two-body potential. 

In tight-binding models the interaction is often written 
in second quantized form as 



H = 



2 ^ 



Vij,i>fbMbi>b f , 



(7) 



where b\ is the usual creation operator of a state <I>i at 
site i, and 



Vii.i' i' 



^(z 1 )^ J (z 2 )^*,(z 1 )^,(z 2 )\^ m (z 1 ,z 2 ,....)\ 2 



V(z% - z 2 )dzidz 2 , 



(8) 



are the four-center integrals. Here, the Laughlin wave 
function present in the total wave function is included in 
which all other particles have been integrated out. I^ml 2 
does not affect the integrals much except that it contains 
a correlation hole ~ (z\ — z 2 ) 2m which screens short dis- 
tances of the two-body interaction. 

The largest of the Vij^y are the interaction between 
particles on just two sites, the direct i = i',j = j' and 
the exchange i = = i'. As described in Appendix A 
the other overlap integrals are small because the local- 
ized wave functions are rapidly oscillation in angle and 
interfere destructively. We therefore extract the largest 
(tight binding) direct and exchange part of the FQHE 
Hamiltonian of Eq. Q 



2 ^ 



(9) 



Here, rij = b\bi is the number operator on site i. 

For long range interactions such that the two-body in- 
teraction can be treated as a constant, V ~ Vb, only 
the direct interaction survives, i.e. Vij.i'y = VoSijiSjj' , 
and therefore H — Hq. This is the mean field limit with 
energy E = VqM b v b - 



3 



For finite range interactions terms will remain in the 
Hamiltonian that connect different orbitals. These will 
be denoted 

T = H-H . (10) 

T contains hopping terms between nearest neighbors 
as the Hubbard Hamiltonian but also contains hopping 
terms to any neighbor state n. As the wave functions are 
rapidly oscillating as described in Appendix A and not in 
phase in T, their overlap is generally smaller as those in 
Hq. We shall in section V discuss the effect of the hop- 
ping terms in T and possible Mott-Hubbard transitions 
if T is strong. 

As the direct and exchange potentials depend only on 
the distance n — j—i between the two interacting particle 
sites, we can rewrite the couplings in Hq as 

U n = -{Vij^ + Vijji) . (11) 

Note that Um b -u — U n . In the following this double 
counting is avoided by truncating n < Mb/2 and multi- 
plying by a factor of two. 

In the following only U n will be important whereas 
the details of the localized wave function and the two- 
body interaction are contained in U n . We shall only as- 
sume that the two-body interactions are repulsive and 
decreasing with range as is the case for Coulomb repul- 
sion between electrons. As result the resulting couplings 
U n will also be repulsive and decreasing functions of the 
neighbor distance, i.e. U n > U n+ \ > 0. 

Fortunately, as will be shown in the following section, 
the Hamiltonian Hq can be solved for rather general cou- 
plings U n and a number of qualitative results can be 
drawn such as incompressible states at fractional fill- 
ings. Whether Hq is a good approximation for the FQHE 
Hamiltonian H as argued above will depend on the effect 
of T = H — Hq on the ground state as will discussed in 
section V. 



III. THE GROUND STATE OF THE EBHM IN 
THE STRONGLY INTERACTING LIMIT 

Eqs. © and l(TT)> result in a ID Bose-Hubbard Hamil- 
tonian in the strongly interacting limit extended with 
far-neighbor interactions U n but without hopping terms 

M B -l M B /2 

H = ^ X! U n n t [n l+n - 5 nfi ] . (12) 

i=0 n>0 

The EBHM also describes optical lattices and reduces to 
Hq for strong lattices in the strongly interacting limit 
near Feshbach resonances 0, ^| . 

The Bose-Hubbard model has only on-site interactions, 
i.e. only Uq is non-zero. The Ising [(| spin model (<Ji = 
rii — 1/2) has only nearest neighbor interaction, i.e. only 
U\ is non-zero. A model including next-nearest neighbor 



interactions, i.e. U2 # was solved by Stephenson |8| and 
Majumdar & Gosh |!|. These statistical models mostly 
investigate phase transitions and critical exponents in one 
or more dimensions for fixed chemical potential. We are 
here interested in the ground state (zero temperature) in 
one dimension in the thermodynamic limit N, Mb — > 00 
for constant density or filling fraction vb — N/Mb- 

We shall investigate a general class of EBHM with 
strong repulsive interactions, i.e. Hamiltonians Hq with 
repulsive U n > 0, n — 0, 1,2,3, etc., that decrease with 
distance such that U n > U n +±, i.e. dU n /dn < 0. Eventu- 
ally we may expect that U n becomes negligible at large n 
as compared to other effects from impurities, finite tem- 
perature or the experimental resolution. In addition, we 
shall assume convexity [Io| 

U n+1 + U n - 1 -2U n >0, (13) 

i.e. d 2 U n /dn 2 > 0. 

The ground state is then organized as to minimize its 
energy. This is achieved by distributing the particles and 
also vacancies as "evenly" as possible avoiding "clump- 
ing". We shall argue that the particles sites (and va- 
cancies) are ordered with a periodicity or sequences of 
replicas depending on the filling fraction (see also [l(j). 

A. Fraction of n'th neighbors, F n 

As an example, consider the filling fraction 5/9. The 
ground state is build up of a series of Mb/9 identical 
sequences (replicas) repeated after each other of the fol- 
lowing type 

I ••0»0»0»0 I ••0»0»0»0 I I/g =5/9 

Here: »=site with a particle and o=hole or vacancy. The 
delimiter | separates replicas. Interchanging any • with a 
o increases the energy since U n > U n +i, and it is there- 
fore the ground state. By translating the whole lattice 
the nine degenerate ground states are obtained. However, 
they any not connected by any two-body interactions. 

This principle of "evenly ordering" generalizes to any 
rational number for the filling fraction N/Mb = q/n. 
i.e. q particles in n states. If q < n there will be (n — 
q) vacancies. The minimum energy state is when the q 
particles (or n — q vacancies) are as distributed over the 
n sites avoiding clumping. 

In the thermodynamic limit N, Mb — > 00 it is now 
a matter of combinatorics to add up the number of n'th 
neighbors. The ground state energy is necessarily a linear 
function of the couplings U n 

Mb/2 

E{N/M b ) - M b Fn(v B )U n . (14) 

The coefficients F n {ys) are functions of the filling frac- 
tion. They are the fractions of sites occupied with an 
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FIG. 1: The fraction of occupied sites with a neighbor n sites 
away, F n (vB), for n = 1,2, 3, 4, 5. 
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FIG. 2: As Fig. 1 but for n = 6, 7, 8, 9. 



occupation also n sites to one side (not to both sides in or- 
der to avoid double counting). In other words MbF^vb) 
is the number of sites occupied which also has a particle 
displaced n sites to one side (only to one side when the n 
and Mb — n double counting is avoided). This fraction is 
necessarily smaller than the filling fraction, F u (vb) < vb- 
The ratio F n {yB)jvB is the conditional probability for an 
occupied site to have a neighbor occupation n sites to one 
side. 



TABLE I: The coefficients in F n = a n v B + P, 





f3„ 


range of vb for q = 0, 1, [n/2] 








[ q/(n-l) , (q+l)/(n+l) ] 


1+n 


-q 


[ q/(n+l) , q/n ] 


1 - n 


q 


[ q/n , q/(n-l) ] 






range of vb for q = n — 1, n — 2, [n/2] 


1+n 


-q 


[ (q-l)/(n-l) , q/n ] 


1 - n 


q 


[ q/n , (q+l)/(n+l) ] 


2 


-i 


[ q/(n+l) , (q+l)/(n+l) ] 



It will be shown below that the functions F n are piece- 
wise linear functions of vb (see also Figs. 1+2) 

F n = a n v B + p n . (15) 

The prefactors a n and (3 n are discontinuous at the frac- 
tions but constant between. They are in fact simple 
integers of changing signs at the fractions vb = q/n, 
q = 1,2, ....,71- 1. 

Before deriving their general dependence on vb ana- 
lytically a few examples illustrates the general idea. At 
half filling the state is 

....•o»o»o»o»o»o»o»o Vb = 1/2 

This is the charge density wave phase of crystalline order 
at half integer density. This incompressible Mott insula- 
tor state also undergoes a transition to a supcrfluid at a 
certain hopping strength 0, 0, 0] . 
Its energy is from Eq. (I14|> 

E(1/2) = N ( 16 ) 

n=2,4,6,... 

i.e., F n = for odd n and F n = 1/2 for even n. 

Let us next study F\ and F2 with increasing filling. 
The fraction of occupied neighbor states, JPi, is zero up 
to half filling, j/b^1/2. Above it increases linearly as 
Fi = 2vb — 1 reaching unity at full filling vb — 1- The 
fraction of next-nearest neighbors is zero for ^s^il/3 and 
then increase linearly as F2 = ?>vb — 1 up to half filling, 
where after it decreases linearly as F2 = 1 — i>b up to 
vb = 2/3. 

...|..o|..o|..o|..o|..o|..o| „ B =2/3 

For filling above vb > 2/3 the number of next-nearest 
neighbors again increase as F2 = 2vb — 1 just as Fi. 

The ordering principle implied by the condition 
£7 n+ x < U n allows us to calculate the probability co- 
efficients F n generally. Let us start at low filling frac- 
tion. For vb < + 1) there are no n'th neighbors 
and F n — 0. The coefficient now increase linearly up to 
vb = where every particle has an n'th neighbor and 
therefore F n — vb — 1/n. By increasing the filling frac- 
tion further the number of n'th neighbors decrease lin- 
early and in fact vanish again at vb = l/(n — 1). F n now 
remains zero from vb — ~ 1) up to vb = 2/(n + 1) 
where the zig-zag behavior repeats. At vb = 2/n every 
particle again has a neighbor n sites away (and also one 
in between at n/2 for even n and at (n± l)/2 for odd n). 
Therefore F n (vB) = vb = 2/n 

This zig-zag behavior of F n in fact repeats itself around 
each vb = q/n, q = 1,2, ...,n — 1 (see table 1 and Ap- 
pendix B) with the same a n and (3 n . Above vb > 1/2, 
however, with F n = 2vb — 1 between the zig-zag's instead 
of zero. 

Because F n is piecewise linear between these frac- 
tions the coefficients can be calculated directly using that 
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F n {q/n) = q/n for u B < 1/2 and i/g > 1/2 respectively. 
The coefficients are (see also Table 1) 



C. Chemical potential gaps and superfluid states 



On 



= l±n 
= TQ 



(17) 



in the intervals [q/n;q/(n± 1)]. Between these intervals 
v B G [q/(n ~ 1)) (q + l)/( rl + 1)] they both vanish for 
vb < 1/2. 

For 1/2 < vb < 1 we find the very same coefficients 
just below and above g/n as in Eq. (|17J) . However, the 
range is different (see Table 1) and between these inter- 
vals F n — 2vb — 1, i.e. a n — 2, (3 n = —1. The coefficients 
are also shown in Figs. 1+2 for n = 1,2, 9. 

For larger filling fractions, vb > 1, the coefficients can 
be related as 



F n (v B ) = F n (v B - 1) + (2v B - I) 



(18) 



i.e., the zig-zag behavior repeats for all F n at vb = 1 + 
q/n. Note that Fo vanishes for vb < 1 but increases 
linearly as Fq = vb — 1 with filling fraction 1 < vb < 2. 



B. Particle- hole symmetry 

The coefficients -F n and resulting energies can be re- 
lated at vb and 1 — v B by exploiting particle-hole sym- 
metry. Since M B F n (y B ) is the number of particles with 
a particle n sites away, Mb{vb — F n (vB)) is the num- 
ber of particles with a hole n sites away. The latter is 
by particle-hole symmetry also equal to the number of 
holes with a particle n sites away, which again is equal 
to M B (1 — vb — F n (l - u B )). Here M B (1 - v B ) is the 
number of holes and MbF„(1 — vb)) the number of holes 
with a hole n sites away. Equating the two numbers for 
the number of particle-hole neighbors n sites away gives 

F n (v B ) = F n (l-v B )+ 2v B - 1. (19) 

By insertion in Eq. I|14f) we obtain the relation 

E(v B ) = E(l - v B ) + (2v B - l)E(l) (20) 

Here E(l) — J2 n >i Un is the energy for N = M B - 

A similar relation has been derived for the ground state 
energy of the FQHE also exploiting particle- hole symme- 
try m 



E(v) = E(l -v)+ {2v - 1)E(1) . 



(21) 



Although the two relations appear identical, they are not 
because v and vb are different. One obvious cause for 
this difference is the neglect of T in the EBHM energy 
of Eq. ((201). For N > M B the tight-binding approxi- 
mation is poor as many states overlap. Furthermore, the 
EBHM energy does not take the Laughlin wave function 
into account when the energy is calculated. This be- 
comes increasingly important for increasing filling frac- 
tion AT — » M because the Laughlin wave function con- 
tains N powers of z whereas the EBHM wave function 
only has Mb = M — N powers of z. For these reasons 
the EBHM ground state energies is not applicable to the 
FQHE for large filling fractions. 



The chemical potential at zero temperature is now 



(22) 



M B 



The chemical potential is therefore also constant between 
the fractions and exhibit plateaus as shown in Figs. 3+4. 
At the fractions the chemical potential is discontinuous 
and therefore an insulator - but a conductor between the 
fractions. A chemical potential gap appears in the spec- 
trum at each fraction q/n, 



A(i(q/n) = n(q/n) + - ^q/n)-, 



(23) 



where ± refers to the right/left side of the incompress- 
ible fraction q/n. The gap can be calculated from the 
expression (|17f) for a n . For 1/2 < i/g < 1 we find 



Au(q/n) = n(J7„_i - 2U n + U n+1 ) 
dn 2 



(24) 



Note that the convexity condition of Eq. (|13fl insures 
that the gap is positive and the system therefore stable. 
For vb < 1/2 wc obtain 

Au{q/n) = (n-2)U n - 1 -2U n + {n + 2)U n+1 
'd 2 U n \ JdXJ n 
dn 2 J \ dn 



(25) 



Again the convexity condition and that dU n / dn < in- 
sures that the gap is positive. 

If q and n have a common divisor as, e.g. 4/6=2/3, 
then the gaps add up. The convexity condition may 
therefore be unnecessarily strict for such fractions. 

If the couplings decrease as a power law U n oc n^ 1 , 
then Afi(q/n) = 7(7 + l)U n /n, and the gaps decrease 
faster than 1/n. Numerical calculations of a finite num- 
ber of electrons on a sphere find FQHE gaps that de- 
crease slightly faster than 1/n |2j|. The couplings U n oc 
log(M B /n) leads to gaps scaling as ~ 1/n. The Coulom- 
bic interactions and wave functions localized in z would 
lead to U n ~ e 2 /l(n + const.) as discussed above. The 
interaction strength may, however, be weakened by layer 
thickness and effects of disorder |2||. The corresponding 
gaps would thus scale as ~ e/ln 2 at large n. 

In the FQHE problem the repulsion between electrons 
is canceled on average by a background field from the 
positively charged solid. This mean field is, however, 
inert and can be ignored in calculated the pressure and 
chemical potential of the quantum fluid. Likewise, the 
optical lattice acts as an inert background of the atomic 
gas systems. 



D. Pressure and incompressible states 

The pressure P = -(dE/dV) = (fiN - E)/V at zero 
temperature can be calculated by noting that the volume 
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FIG. 3: The energy, chemical potential and pressure vs. filling 
fraction vb- A short range two-body interaction is chosen 
with neighbor couplings U n = 3/(nn) 2 . 



FIG. 4: As Fig. 3 but with neighbor couplings U n = 
log(A^B /n) /2Mb for n < Mb = 20, which simulates a longer 
range two-body interaction. 



is proportional to the number of states Mb [2lJ . We shall 
simply set the volume per state to unity whereby 

f dE \ E v-^ 

The pressure is therefore constant, i.e. compressible, be- 
tween the fractions. If the couplings U n are convex we 
observe that the pressure jumps to a higher pressure and 
the system is therefore an incompressible quantum fluid 
at those fractions as shown in Figs. 3+4. If the couplings 
U n are not convex at n, the pressures may decrease at 
vb = q/n+ which would require a Maxwell construction. 
The pressure discontinuities 

AP(q/n) = P(q/n+) - P(q/n_) , (27) 

can likewise be calculated by inserting the (3 n of Eq. i|17|) 
in the pressure. For vb < 1/2 we find 

AP{q/n) = q(U n -i - 2U n + U n+1 ) 

assuming that U n is a smooth function of n. For 1/2 < 
vb < 1 we find 

AP(q/n) = (q-l)U n - 1 -2qU n + (q+l)U n+1 (29) 



E. Gap dependence on two-body range 

In Figs. 3+4 the energy, pressure and chemical poten- 
tial are shown as function of filling fraction vb for two 
representative set of couplings. 

The first assumes a short range (delta function) two- 
body interaction such that the couplings are U n = 



3/ (7m) 2 as discussed above. These couplings decrease 
rapidly with neighbor distance n and therefore the chemi- 
cal potential gaps and pressure discontinuities are largest 
at the fractions q/n with small n, e.g. q/n = 1/2, 1/3, 
2/3, 1/4, 3/4, etc. 

In Fig. 4 a longer range two-body interaction is as- 
sumed such that the neighbor coupling decrease more 
slowly with neighbor distance n. The specific choice is 
U n = \og{M B /n)/2M B for n < M B = 20 and van- 
ishing coupling for larger n. A logarithmic coupling 
is obtained by averaging a two-body Coulomb poten- 
tial over radial variables preserving a phase difference 
4>i — 4>2 — 2ttti/Mb- The rather slow decrease of U n 
with n has the effect that the gaps and pressure discon- 
tinuities only decrease as *~ l/n. The long range inter- 
action approaches the mean field limit V = Vq where 
E/M B = v%/2 = P and fi = v B - 



IV. INCOMPRESSIBLE STATES AT 
FRACTIONAL NUMBERS 

The incompressible states found for the system of 
bosons above can as described in section II be translated 
to the FQHE states for fermions if we assume that the 
hopping term T in the full Hamiltonian is small, or if it 
does not destroy the incompressible state by making a 
Mott transition. 

The incompressible states in the EBHM appear at frac- 
tions vb = q/n, where n and q are integers. However, 
when applied to the FQHE, it is important first to sepa- 
rate out the Laughling wavefunction with maximal m as 
mentioned in section II. For this reason values vb < 1/2 
do not come into play in the FQHE because the pre- 
ferred ground state wavefunction then includes a Laugh- 
lin wavefunction with a to- value increased by 2. 

According to Eq. the Bose fractions correspond to 
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TABLE II: The leading incompressible filling fractions for 
small n and q for m = 1 (see text). Observed FQHE fractions 
are given in bold face. 



n 


q 


v = vb/(1 + mva) — q/(n + mq) 


1 


1,2,3,4,5,6,... 


1/2, 2/3, 3/4, 4/5, 5/6, 6/7,... 


2 


1,3,5,7,... 


1/3, 3/5, 5/7, 7/9,... 


3 


2,4,5,... 


2/5, 4/7, (5/8),... 


4 


3,5,7,... 


3/7, 5/9, 7/11),... 


5 


3,4,6,7,.. 


3/8, 4/9, 6/11, 7/12, .. 


6 


5,7,11,... 


5/11, 7/13, 11/17,... 


7 


4,5,6,8,... 


4/11, (5/12), 6/13, 8/15, ... 


8 


5,7,9,11,... 


5/13, 7/15, 9/17, 11/19,... 


9 


5,7,8,10,11,.. 


(5/14), 7/16, 8/17, 10/19, 11/20,... 


10 


7,9,11,13,.. 


(7/17), 9/19, 11/21, 13/23,... 


11 


6,7,... 


6/17, 7/18,... 



should be employed. For example, for 1/5 < v < 1/3 
(corresponding to 1/4 < vb < 1/2 for m = 1) the Laugh- 
lin wave function with to = 3 should be employed instead 
with a larger vb such that v — vb/(3vb + 1). 

For to = 3 the commensurable fillings vb = n, 
n = 2, 3, 4, 5, therefore lead to the fractions 2/7, 3/10, 
4/13, 5/16, etc. These first three are clearly identified 
in the FQHE spectrum |3j. The next leading fractions 
v B = 3/2, 5/2, 7/2, etc. gives v = 2/9, 5/17, 7/23, etc., 
respectively. The two first are also observed in the FQHE 
spectrum. The large vb become gradually more difficult 
to observe because the corresponding FQHE fraction ap- 
proaches 1/3 and is therefore quenched by the Laughlin 
fraction v — 1/3. 



FQHE filling fractions 



VB 



q 



l 



mvB 



n + mq 



(30) 



where to = 1,3,5,... is the Laughlin number. Most of 
the leading fractions are shown in Table II and will now 
be compared to the FQHE fractions found in FQHE ex- 
periments @,0|. 



A. The Laughlin and Jain sequences 

The incompressible states at vb = 1/2 gives the 
Laughlin fractions 1, 1/3, 1/5, 1/7, etc. for m = 1, 3, 5, 7, 
... respectively. Next to these the most important set of 
fractions is the Jain sequence 



2n± 1 



(31) 



These appear for m = 1. Inserting the EBHM fractions 
vb = 1 — 1/n, n = 1,2, 3, ... in Eq. gives the positive 
(+) part of the principal FQHE fractions, i.e. 1/3, 2/5, 
3/7, etc. There are analogous incompressible states at 
fillings vb = 1 + 1/n. Inserting these in Eq. Q they 
lead to the other (-) half of the principal fractions of Eq. 
HIJ, i.e. 2/3, 3/5, 4/7, etc. 

The Jain fractions will be leading for two reasons. 
Firstly, the pressure discontinuity of Eq. (|28[1 increase 
with q and q = n — 1 is maximal for these fractions for a 
given n. Secondly, in the region above vb = 1 — l/(n— 1) , 
the fraction vb — 1 — 1/n will be the incompressible frac- 
tion with the smallest n. As the pressure discontinuity 
of Eq. I|28[) also decreases with n it will therefore be the 
largest in this region. 



B. 



v < 1/3 



It was discussed in section II that for fractions v < 
1/to, to = 1,3,5, .. the Laughlin wave function with to 



C. 1/3<i/<1/2 

Secondary incompressible states also appear in the 
EBHM at the fractions q/n, q = n — 2,n — 3,..., etc. 
As their gap size decrease with n their importance are 
ordered likewise (see Table II). The leading secondary 
fractions (not including the Laughlin and Jain series) are 
therefore: v B = 3/5, 4/7, 5/7, 5/8, 5/9, 7/9, 7/10, 6/11, 
etc. These EBHM fractions correspond to the FQHE 
fractions v = 3/8, 4/11, 5/12, 5/13, 5/14, 7/16, 7/17, 
6/17, etc. (see also table II). Remarkably, four of these 
fractions are clearly identified in recent high resolution 
FQHE experiments |3|. The fraction at 5/12 is compat- 
ible with the structure reported between the principle 
fractions v = 2/5 and 3/7. The fractions 5/14 and 7/17 
are compatible with structures between 4/11 and 6/17. 
The remaining fraction 7/16 is squeezed between the two 
Jain fractions 3/7 and 4/9 and thus difficult to observe. 

Incompressible states with fractions q/n with larger n 
will be weaker. Except for the Jain fractions, they will 
also be placed between stronger incompressible states, i.e. 
fractions with smaller n. Therefore they will be increas- 
ingly difficult to identify experimentally. The incom- 
pressibilities and chemical potential gaps become smaller 
and eventually disappear with increasing n. Effects of 
impurities, finite temperature and experimental resolu- 
tion therefore exclude experimental observation of the 
incompressible states at large n. 



D. 



v = 1/2 



The EBHM is apparently successful in producing all 
the incompressible states at the principle and secondary 
FQHE fractions including those found in recent experi- 
ments. However, it clearly fails for m = 1 at half filling 
or vb = 1 where it also predicts an incompressible state 
in contradiction with FQHE experiments and theoretical 
calculations [24| . In the Bose-Hubbard model the state at 
N = Mb becomes compressible for a sufficiently strong 
hopping term. It will be discussed in the next section 
that the hopping term T may cause a transition from 
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a Mott insulator to a superconductor, i.e. changes the 
state at v — 1/2 to a compressible one. 



E. v>\j2 

For vq > 1 the EBHM predicts a number of incom- 
pressible states. The leading ones are shown in Table II. 
Experimentally, they are increasingly difficult to observe 
as v — ► 1 due to limited resolution. The set of states at 
v B = 2- l/n = 5/3, 7/4, 9/5, etc., correspond to FQHE 
states v = 5/8, 7/11, 9/14, respectively. The 7/11 is 
clearly observed whereas a structure is found at 5/8 and 
possibly also at 9/14. 

However, as mentioned above the ground state energies 
of the EBHM become less reliable with increasing filling 
fraction as a model for the FQHE. When Vb > 1 it is 
therefore a better approximation to use the particle-hole 
symmetry of the FQHE, Eq. by inserting the EBHM 
energies for vb < 1 as an approximation for the FQHE 
ground state energies. It follows that the incompressible 
states for v also exist for 1 — v. 



V. MOTT-HUBBARD TRANSITIONS 

The term T contains a number of hopping terms be- 
tween n'th neighbors. It is therefore instructive to com- 
pare to the one dimensional Bose-Hubbard model (Eq. 

which contains nearest neighbor hopping t\ but only 
on-site interactions, i.e. only Uq is non- vanishing. In the 
strongly interacting limit (T <C Hq) it is a Mott insulator 
for integer filling factors [ill ll2T | as also found above for 
the EBHM. Increasing the hopping term a Mott transi- 
tion takes place at t\ ~ Uq/\.92 in one dimension and 
the system becomes a superfhhd. 

It is known from studies of the FQHE state at half 
filling (or N/Mb = 1) 2J| that the state is compress- 
ible. We should therefore expect that by including T 
the incompressible state at N = Mb undergoes a Mott 
transition. 

Adding terms U n ^ but keeping T — 0, creates 
new Mott insulator states at the corresponding fractional 
filling factors as described above. A crucial question is 
therefore whether these Mott insulator states survive as 
the hopping term T is gradually turned on. Studies of 
the Bose-Hubbard model including next neighbor inter- 
actions show that the Mott insulator state at half filling 
or charge density wave disappear before that at unity fill- 
ing as the nearest neighbor hopping t\ is increased. T is, 
however, very different from t%. 

A closer look at the ground states for fractional and 
integer filling fractions reveals a difference in the overlap 
between the degenerate ground states. The ground state 
for N = Mb — 1 



is Ms-fold degenerate and the hopping term T couples 
them all. As a consequence the energy can change con- 
siderably. It is the far-neighbor hopping terms in T that 
can be responsible for "destabilizing" the integer filling 
fraction states. Only the nearest neighbor hopping t\ is 
present in the Bose-Hubbard model. 

In comparison the states near fractional filling, e.g. 
N = 2Mb/3 - 1 

....|..o|..o|..o|.o.o|..o|..o|..o|..o| 

this state only couple to two of the N degenerate ground 
states, namely those where the extra hole hops one site 
to one side, resulting in a modest energy corrections. 

Another difference appears when comparing the cor- 
rections to the ground state energies just above and be- 
low N/M B =2/3 with those at N = M B - In the former 
case the correction from including hopping terms from T 
is the same for N — 2Mb/3 ± 1. However, the leading 
hopping term correction for N = Mb ± 1 differs because 
it vanishes when one particle is added to commensurate 
filling. 

It should be noted that only the odd commensurate 
states vb = 1,3,5,... should undergo a Mott transition 
and only for m = 1 in order to fully agree with the FQHE. 
It has be argued that the far hopping terms are impor- 
tant for commensurate states but at present the EBHM 
cannot explain why only the odd commensurate states for 
to = 1 vanish for the FQHE. It may be related to the dif- 
ference in energy between the energies of the EBHM and 
FQHE, Eqs. (19) and (20) respectively. The Laughlin 
wavefunction may affect energies of the FQHE especially 
at large filling vb and to. 

The EBHM is at most an approximation to the FQHE 
since the localized states are only a good approxima- 
tion when the hopping terms are small. Furthermore, 
a bosonization procedure was applied without including 
the Laughlin wave function, when energies were calcu- 
lated within the EBHM. These approximations should 
be investigated further. 

The effects of far-neighbor hopping is a complicated 
yet important problem. Experimentally, it may be stud- 
ied by quantum gases in optical lattices where the neigh- 
bor couplings and hopping can be controlled by tuning 
the lattice spacing and height, interaction strengths and 
particle density. Both the neighbor couplings U n and 
the hopping coefficients t n decrease exponentially with 
the tunneling factor ~ enp(—2n^/Vo / Er) , where Vb is 
the optical lattice amplitude and Er the recoil energy 
[l2j |. In addition, however, the couplings U n also scale 
with the scattering length which allows us to enhance the 
couplings near a Feshbach resonance and thus reach the 
strongly interacting limit. If absorption images can be 
obtained from such a one dimensional optical lattice then 
Bragg peaks should be observed displaying fractional re- 
ciprocal lattice vectors due to the periodicity at fractional 
filling. 
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VI. SUMMARY AND OUTLOOK 

The Bosc-Hubbard model extended by including near- 
est and far neighbor interactions was related to the frac- 
tional quantum Hall effect problem. The ground state 
is calculated for the extended Bose-Hubbard model in 
the strongly repulsive limit (weak hopping) . Incompress- 
ible Mott insulator states were found at filling fractions 
compatible with the FQHE filling fractions of the lowest 
Landau level. The incompressible states with the largest 
gaps are the Laughling and Jain fractions. Secondary in- 
compressible states occur for filling fractions 3/8, 4/11, 
5/13, 6/17, 2/7, 3/10, 4/13, 5/16, 2/9, 5/17. Such frac- 
tions have been observed in recent FQHE experiments. 
Weaker states at 9/14, 5/12, 5/14, etc. are also pre- 
dicted and indications of structures at the fractions are 
reported in recent high resolution FQHE experiments. 
Furthermore, incompressible states at fractions with even 
larger denominators are also predicted; the next candi- 
dates would be 7/16, 7/17, 7/23, etc. However, these 
states are increasingly difficult to observe and is a chal- 
lenge to find experimentally. 

The EBHM provides an interesting spectrum of in- 
compressible state at fractional fillings. Whether these 
states survive or undergoes a Mott insulator transition 
to a superfluid as the hopping terms are included is 
an important question. It has been indicated that far 
hopping terms are important for commensurate states 
vb = 1, 2, 3, ... but it is not clear why only the odd ones 
should undergo a Mott transition. The EBHM is at most 
an approximation to the FQHE since the localized states 
are only a good approximation when the hopping terms 
are small. Furthermore, a bosonization procedure was 
applied without including the Laughlin wave function, 



when energies were calculated within the EBHM. These 
approximations should be investigated further. 

In spite of its limitations as a model for the FQHE, 
the EBHM does provide a ground state wave function in 
the strongly interacting limit from which incompressible 
fractions, gaps and incompressibilities can be calculated. 
Corrections from hopping terms may be estimated per- 
turbatively using this wave function as a starting point. 
Therefore, the EBHM may provide a useful basis for also 
calculating a number of transport coefficients such as the 
longitudinal and transverse Hall resistances as well as 
their low temperature dependence. 

The Mott transitions at fractional fillings between in- 
sulator and superfluid phases may generally be studied 
by experiments with quantum gases in optical lattices. 
By lowering the lattice heights the overlap between near- 
est and next-nearest neighbors increase and therefore also 
U\, U^,..., and the corresponding hopping terms j^. In 
optical lattices the neighbor couplings and hopping can 
be controlled by tuning the lattice spacing and height 
as well as the interaction strengths between atoms and 
the particle density. Furthermore, Bose and Fermi atoms 
can be mixed resulting in a multitude of Mott, c harg e and 
spin density wave, superfluid and other phases |28l l29| . 

In high temperature superconductivity (HTc) or- 
dered charge density waves are found in the spin glass 
phase. These are ordered ID stripe and 2D checker- 
board phases with magic hole doping fractions at x = 
1/16,3/32,1/8,3/16 |30|. These phases and their com- 
petition with HTc are believed to be described by (ex- 
tended) 2D Hubbard Hamiltonians and therefore such 
models are actively investigated. The ordering, Mott in- 
sulator phases and transitions in lattices with fermions 
and/or bosons is an important topic in many fields. 
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APPENDIX A: LOCALIZED STATES ON A 
CIRCLE 

To motivate the basis of localized wave functions in 
Eq. (JSJ we discuss the model of Dyakonov [22| , where the 
FQHE problem is approximated by putting the states on 
a circle. 

The idea behind the circle approximation is to replace 
the coordinates z = re 1 ^ by their phase: e 1 ^, by approx- 
imating the radius by its average value: (\z\ 2 ) = (r 2 ) — 
k + 1, in units of the magnetic length. The states now 
become: ipk ~ exp(ikcf>), k — 0,1,2, Mb — 1- The 
Laughlin wave function is then 

* m -Aj](e' ;0, -e^) m , (Al) 

with m — 1,3,5,... and normalization A 2 — 
(m!/2 7 r) JV /(mAT)!. 

The Wannier-type localized functions analogous to (J3J) 
are simply 

Mb-1 

1 1 — e iMB ^ 

where j = 0,1, ...,Mb — 1. These wave functions are 
localized in angle <fi around 2-Kj/M B - They oscillate 
rapidly with frequency 1/Mb and amplitude dropping off 
slowly as j . They are the same function (a phase times 
s'm(MBx)/ sin(a;) with x — <p/2 — irj /Mb) of the shifted 
angle (4> — 2irj /Mb)- The localization sites are relative 
to each other and can all be shifted by an arbitrary an- 
gle. A many-body wave function build of such localized 
states would therefore rather resemble a strongly corre- 
lated liquid than a crystal. 

On the circle the Vij^if are simply 

\* m (<k,(h, -)\ 2 V(0i - <h)*h.dfa , (A3) 

where now the two-body potential V(4>i — 4>j) depends 
on the angular distance instead of \zi — Zj\ as in Eq. JSJ. 



The relative simple localized wave functions allows us 
to study the overlap integrals in Vij^/ji. The localized 
wave functions are rapidly oscillations and generally out 
of phase except for the coupling between two sites. As 
result the direct coupling Vij^j and the exchange Vijji are 
the largest. We therefore expect that the tight binding 
approximation of Eq. @ , which includes only the direct 
and exchange interactions, is a good approximation. This 
is especially the case for long range interactions which 
approaches the mean field limit. In this limit Vij,i'j' = 
$i,i'8j,j> fV((f>)d(f>. 

For interactions of range shorter than £l/M B , the two- 
body potentials can be approximated by a delta-function 
interaction: V(<f>) ~ 27rVo<5(0). The resulting U n for the 
delta function interaction are 

U n =2TTV \<5> t {<t>)\ 2 \<S> l+n {<j))\ 2 d<t> . (A4) 
Jo 

Here, we have ignored the Laughlin wavefunction but 
will discuss the effect of the correlation hole below. Note 
that generally U n — Mb J V(4>)d4> and for the delta- 
function interaction J2 n U n = 2ttMbVq. Also U n = 

Um b -u- 

Inserting the wave functions of Eq. I|A2(I in (|A4|) yields 
U = (2/3)M B V and U n = M B V [M B sm{irn/M B )r 2 
for n > 1. In the thermodynamic limit Mb — > oo the 
couplings scale as U n = MbVqi '(ra) 2 . For longer range 
interaction the couplings U n decrease with n at a slower 
rate. For a Coulomb-type potential V — e 2 /\4>i — 4>z\ we 
obtain U n ~ Mge 2 /{n + const.), where the constant is 
of order unity and reflects that U n is finite due to the 
correlation hole in the Laughlin wave function. 

The delta-function approximation for the short range 
interaction should only be employed for the bosonic part 
of the wave function. As a consequence the direct and 
exchange interaction energies are the same and leads to 
a factor two. The delta function interaction should not 
be employed for the fermionic part of the wave function. 
The exchange part would then cancel the direct part and 
lead to vanishing interaction energy as it should for an 
antisymmetric wave function. 

For a longer range potential the exchange and direct 
parts do not cancel. Also the couplings U n will decrease 
slower with n. Such an example U n oc log(M B /n) is 
shown in Fig. 4. 

Returning to the Coulomb potential in coordinate 
space V = e 2 /\z\ — Z2\ and the localized wave functions 
of Eq. 101, the couplings of Eq. l(TT|) with (JSJ become 
U n ~ (e 2 /( r )) log(Mg/n). These wave functions are, 
however, only localized in their phases. Therefore the 
2-D integration over both radius and angle is different 
than the ID integration over angle on the circle. The 
result is a log(Ms/n.) in stead of 1/n respectively. 

It may be advantageous to try to construct a basis of 
wave functions that are localized both radially and in 
angle 0, i.e. in localized in z = re 1 ^ . For a Coulomb po- 
tential the couplings will also be approximately Coulom- 
bic U n ~ e 2 /l{n + const.). Here, the constant takes into 



11 



TABLE III: The coefficients in F n = a n v B + f3 n for n < 4. 



distance n 


a n 


/3„ 


range of vb 











[0,1] 




1 


-1 


[1.2] 


1 








[0,1/2] 




2 


-1 


[1/2,1] 


2 








[0,1/3] 




3 


-1 


[ 1/3 , 1/2 ] 




-1 


1 


[ I/ 2 , 2/3 ] 




2 


-1 


[2/3,1] 


3 








[0,1/4] 




4 


-2 


[ 1/4 , 1/3 ] 




-2 


1 


[ 1/3 , 1/2 ] 




4 


-2 


[ I/ 2 , 2/3 ] 




-2 


2 


[ 2/3 , 3/4 ] 




2 


-1 


[3/4,1] 


4 








[0,1/5] 




5 


-1 


[ 1/5 , 1/4 ] 




-3 


1 


[ 1/4 , 1/3 ] 




o 


o 


f 1 /3 9/5 1 




5 


-2 


[ 2/5 , 1/2 ] 




-3 


2 


[ 1/2 , 3/5 ] 




2 


-1 


[ 3/5 , 2/3 ] 




5 


-3 


[ 2/3 , 3/4 ] 




-3 


3 


[ 3/4 , 4/5 ] 




2 


-1 


[4/5,1] 



account that the coupling Uq is finite according to defi- 
nitions in Eq. (JHJ). 



APPENDIX B: TABLES OF F„ COEFFICIENTS 



The first coefficients F n = a n VB + (3 n are listed in 
Tables III+IV for filling fractions vg < 1 and n < 7. 



TABLE IV: The coefficients in F n = a n v B +f3 n for n = 5, 6, 7 
in the interval of fractions 1/2 < vb < 1- 



n 


a n 


/3n 


range of vb 


5 


6 


-3 




1/2,3/5] 




-4 


3 




3/5 9/3 1 




2 


-1 




2/3 , 3/4 




6 


-4 




3/4 , 4/5 




-4 


4 




4/5 , 5/5 ] 




2 


-1 




5/6 11 

/ 5 J 


6 


-5 


3 




1/9 4/7 1 




2 


-1 




4/7 , 3/5 




7 


-4 




3/5 9/3 1 




-5 


4 




2/3 , 5/7 




2 


-1 




5/7 , 4/5 




7 


-5 




4/5 , 5/6 




-5 


5 




5/6 , 6/7 ] 




2 


-1 




6/7 ] 1 ] 


7 


8 


-4 




1/2 , 4/7 ] 




« 
~u 


4: 




4/7 , 5/8 ] 




2 


-1 




5/8 , 2/3 ] 




8 


-5 




2/3 , 5/7 ] 




-6 


5 




5/7 , 3/4 ] 




2 


-1 




3/4 , 5/6 ] 




8 


-6 




5/6 , 6/7 ] 




-6 


6 




6/7 , 7/8] 




2 


-1 




7/8 , 1] 



